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Abstract 

We extend previous calculations of the zero temperature superfluid fraction fs (SFF) vs local- 
ization, from the fee lattice to the experimentally realized (for solid ^He) hep and bcc lattices. 
The superfluid velocity is assumed to be a one-body function, and dependent only on the local 
density, taken to be a sum over sites of gaussians of width a. Localization is defined as cr/d, with 
d the nearest-neighbor distance. As expected, for fee and bcc lattices the superfluid density tensor 
is proportional to the unit tensor. To numerical accuracy of three-places (but no more), the hep 
superfluid density tensor is proportional to the unit tensor. This implies that a larger spread in 
data on fs, if measured on pure crystals, is unlikely to be due to crystal orientation. In addition, 
to three decimal places (but no more) the curves of fs vs a/d are the same for both the hep and fee 
cases. An expected value for the localization gives an fs in reasonable agreement with experiment. 
The bcc lattice has a similar curve of fs vs a/d, but is generally smaller because the lattice is more 
dilute. 



PACS numbers: 67.40.W, 67.80.S 



I. INTRODUCTION 



Theories of superfluidity in solids date to 1969, with the work of Andreev and Lifshitz 
on the long-wavelength, low-frequency behavior of such systems,— and to 1970, with the 
work of Leggett on the non-classical rotational inertia (NCRI).^^ In equilibrium, the NCRI 
fraction (NCRIF) is essentially the superfiuid fraction (SFF), denoted /s-— Rotation of the 
walls of an annulus causes the many-body wavefunction to develop a boundary condition in 
the rotating frame wherein a non-integral phase change, associated with backfiow, accrues 
when each particle coordinate is transported around the annulus. An upper bound was 
developed for the NCRIF using a phase-function (f){r} that included particle correlations 
only in a dependence on the local particle density p(r): what in the present work we call a 
one-body phase function. We assume that the particles making up the solid are bosons. 

Although solid ''He takes on the hexagonal close-packed structure, Ref. |^ argued that the 
results would be similar for a face centered cubic structure, and fg was then calculated for a 
density taken to be a sum of gaussians at each lattice site, with lattice constant a fee- This 
gave fs as a function of the localization of the gaussians, taken to vary as exp(— 
Thus fs was given as a function of h/ajcc- An equivalent notation uses the nearest-neighbor 
distance d = afcc/V^ and the rms width a = hj\pl of the Id gaussian exp(— so that 
/s can also be thought of as a function of a/rf, a lattice-independent quantity that is relevant 
when later comparison is made with hep and bcc lattices. 

Recently, Kim and Chan have observed non-classical rotational inertia (NCRI), suggested 
by Leggett in 1970, for solid ^He in Vycor, porous gold, and in bulk.-^^ Given that the 
extrapolated values of the T = NCRIF lie in the same 1-2% range, despite the different 
sample porosities, it seems unlikely that it is a surface effect. Recently NCRI has also been 
observed in solid para-hydrogen (H2),'^ an indication that the effect occurs for all bosons that 
are sufficiently quantum in nature, as measured, for example, by the (dimensionless) de Boer 
parameter A = /j,^/[crLj(me)], where a^j is the characteristic length and e is the characteristic 
energy associated with the Lennard- Jones 6-12 potential V = 4:e\(cTLj/ry^ — (cr^^j/r)^]. 

More recent and better converged calculations based on Refo yielded that, for the ex- 
pected localization ratio from 1976, the superfiuid fraction would be about 2%,«^ in reason- 
able agreement with the observations of Kim and Chan. Therefore in the present work we 
have calculated fs for an hep solid whose density is a sum of gaussians at the lattice sites, 
again within the context of a one-body phase function. To characteristic accuracy of three 
decimal places (but no more), fs for the hep solid is isotropic, and to the same degree of 
accuracy (but no more) the curve of fs vs a/ d has the same form as for the fee solid. On 
the other hand, to the same accuracy the hep superfiuid response is not literally isotropic, 
and not literally equivalent to the corresponding fee lattice. We also obtain for the bcc 
lattice, whose superfiuid density, hke that for the fee lattice, is isotropic. 

These results imply that pure crystals of hep "^He should give, to an accuracy of better 
than a percent, the same fs for all flow directions or, equivalently, for azimuthal flow along a 
fixed experimental cell within which different crystals are grown. On the other hand, if the 
crystal is not perfect, boundary effects may cause the superfiuid fraction to change. This 
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may be the origin of the spread in observed values of the NCRIF.— 

Although the density of real hep ^He is certainly not a sum of gaussians, we expect the 
result that fs is isotropic for the hep lattice is robust. It is not obvious why fs should be as 
very nearly isotropic as it is, although many properties of electronic systems are very similar 
for the fee and hep systems. 



II. METHOD 



Ref|^ showed that a one-body phase function (f) minimizes the flow energy E = 
J dfp{r)v^, where Vg = ;^V0, when the equivalent of the continuity equation is sat- 
isfied: 



h 

m 

For a periodic system we take 

m 



= V ■ J(f) = ■ (p(r) V0(f)) . (1) 



r) = — Wo ■ + ^</)Gexp(zG ■ f), (2) 



so 

Vs{r) =tTo + J^^^Gexp(2G-f) = H ^ G0Gexp(iG ■ r), (3) 

— * 

where Vq is considered to be a known imposed average superfluid velocity and the G's are 
reciprocal lattice vectors for the hep lattice. To find the unknown flow pattern (equivalent 
to finding the unknown phases we use (1) to set up a set of linear equations with vq as 
the source term. From the phase we can calculate the superfluid velocity and the superfluid 
fraction. 

The fourier transforms satisfy 

p{f) = pg exp{iG ■ r), Pg = / dfexp{iG ■ f)p{f), (4) 

G 

where V is the system volume. Explicitly, the convolution theorem gives for the fourier 
components of j(r) = p{r)Vs{r) that 

3g = ^Pg~g'Vg' = ^ Pg-g'Vg' + PgVq. (5) 

G' GVO 

Let 

i^G = — 0G, (6) 

m 

so Vg = GipG- Then the continuity equation implies, for G 7^ 0, that 

= G • jG = ^ Pg-g'G ■ vg' + PgG ■ ^^o = ^ Pg-g'G ■ G'tpG' + PgG ■ vq. (7) 

GVO GVO 
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When the term in vq is brought to the other side of the equation, physically this means that, 
along any G ^ 0, the flow due to the fourier components is opposite the average flow vq. 
This induced flow opposite to Vq also holds for G = 0, since typically < 1. 

The V'g's that determine the flow pattern can be solved for and then substituted into the 
equation for the average flow current 

jo = ■ ^0 = Povo + ^ P-gVg = PoVo + ^ P-gGiPg- (8) 

This completely defines the components of the superfluid density tensor. We obtain its 
component along vq via 

Vq - ps ■Vo = Po + P~G — = Po + P^gVg — —■ (9) 

G^O ^0 G^O ^0 

These equations are valid for any pc- Sometimes we will refer to 

/. = '-^^-^ (10) 
Po 

as the superfluid fraction, although properly we should refer to 

7:=^ (11) 

Po 

as the superfluid fraction tensor. 

We represent a Gaussian in three-space, normalized to unity, by 

p{f) = {^/exp[-C-n (12) 

In Id this has rms width a = b/y/2. As already noted, Refsj^ and 10, for fee lattices plot 
fs as a function of b/ajcc = b/{d\/2) = a/d. In the present work we also consider the hep 
lattice and the bcc lattice. 

Let Po = I/Vq be the number density in terms of the unit cell volume Vq. Then for the 
sum-of-gaussians model the fourier transform of the total density is given by 

Pg = Po exp(-|Gp6V4) cos(G ■ u). (13) 



III. HCP LATTICE 



A. Lattice 

For the hep lattice we take the hep lattice constant in the x-y plane to be given by atcp, 
so the nearest neighbor distance d = atcp- Thus cr/ahcp = cr/d. The volume of a unit cell is 
Vo,hcp ~ ^^hap — V^d^', as can be determined from Vo^hcp = di ■ d2 x a^, with the basis set 
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vectors given below. This is the same as for the fee lattice with the same nearest-neighbor 
distance d = ajcc/V^, since Vqjcc = ct/cc/^ = V^d'^. 

The hep real-space basis is taken to be (omitting the subscript on a) 

ai=a(l,0,0), 02 = ^(1,73,0), 03 = 0^1(0,0,1). (14) 

The lattice sites are specified hjR = idi + ja2 + ka^,, where z, j, k are integers. 
Associated with this is an atomic basis set at iw, where 

The atoms are located at i? ± m, where 2u has length d = a. This choice is not unique; we 

could also have employed u = |(1, -1-;^, \J\)- 
The reciprocal space basis is given by 

6i = — (l,--=,0), 62 = — (0,-=,0), 63 = — (0,0,./-). (16) 

Reciprocal lattice vectors (RLVs) are specified by G = mhi + 77.62 where m,n,p are 

integers. 



B. Results 

In order to obtain numerical results one must truncate the sum over RLVs. To test for 
convergence we computed as the RLV basis set was increased. We employed two methods: 

(1) With 771,71, p all running from — to +A^, we computed fg vs N for each value of 
b/d = \^a/d, and each of the fiow directions considered, which included (1,0,0), (0,1,0), 
(0,0,1), (1,1,1), and (1,-1,0). The results for b/d = 0.35 {a/d = 0.2474) and (100) and (001) 
are given in columns two and three of Table I. The largest value taken for any b/d value was 
N = 8. 

(2) Because the RLV length for the z direction is relatively short, to get a better repre- 
sentation of fiow along z we often included more fourier components in that direction, so 
that p ran from — 2A^ to 2A^. This led to improved convergence rates for fs, denoted by in 
columns four and five of Table I. For this method the largest value taken for any b/d value 
was N = 7. 

Comparing the (most converged) = 6 rows of Table I, we conclude that, although in 
principle fs is a tensor quantity with 2,3. = fg^yy ^ fs,zz, for the hep lattice it may be treated 
as a scalar. Moreover, to the same accuracy, fs for the hep and fee case (columns six and 
seven of Table II) are indistinguishable (actuallxi "we find fs^J^. < f^'^^ < fsTz)i ^^us the 



form of fs vs a/d is the same as given by RefllOl where {b/a)fcc = cr/d. For completeness. 



we present fs vs a /d m. Figure 1, along with data for the bcc lattice, which is discussed in 
the next section. 
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TABLE I: Superfluid fraction fg for b/d = 0.35 (cr/d = 0.2474) and various flow directions. For 
the hep lattice, columns two and three are by method (1) and columns four and five by method 
(2), as discussed in the text. For the fee lattice, columns six and seven are for two different flow 
directions, where a = (lll)/\/3- 





hcp(l) 


hcp(l) 


hcp(2) 


hcp(2) 


fee 


fee 


N 


fs,xx 


fs,zz 


f 


f 


fs,xx 


fs,aa 


2 


0.91992224 


0.923652853 


0.919241677 


0.919773502 


0.91950571 


0.919721444 


3 


0.91910979 


0.920556418 


0.919038558 


0.919612587 


0.91926078 


0.919278107 


4 


0.91903729 


0.919752869 


0.919024750 


0.919607992 


0.91924305 


0.919244375 


5 


0.91902573 


0.919630510 


0.919023744 


0.919607848 


0.91924173 


0.919241833 


6 


0.91902397 


0.919611569 


0.919023668 


0.919607842 


0.91924163 


0.919241638 




FIG. 1: Superfluid fraction vs a/d for the fee, hep, and bee lattices. 



To show the extent that fs is not isotropic, we present Table II for fg^zz — fs,xx, which is 
typically down from fg^^z by three orders of magnitude. The smaller the value of cr/d, the 
less converged the results. 

For converged values of fs, we have verified that, for flow along, e.g., (Ill), the current 
can be obtained by application of jo = ■ "^^o and the tensor components ps,xx and ps,zz 
obtained from flow along (100) and (001). 

We have also employed (8) to obtain the three components of the current flow, for various 
directions of average superflow. These results are consistent with our studies indicating that 
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TABLE II: Anisotropy fs,zz — fs,xx for hep lattice for various values of b/d. Also given is the largest 
value of A'^ needed to obtain convergence (when possible) , using method (2) (see text) . 



b/d 


0.70 


0.65 


0.60 


0.55 


0.50 


0.45 


0.40 


0.35 


0.30 


0.25 


0.20 


{fs,zz fs,xx) ^ 10 


0.001 


0.004 


0.019 


0.062 


0.159 


0.321 


0.495 


0.583 


0.517 


0.330 


0.136 


N 


2 


3 


3 


3 


3 


3 


4 


6 


7 


7 


7 



the superfiuid density tensor is nearly isotropic. Moreover, it is isotropic in the x-y plane 
when the calculations are performed to convergence, just as the full superfiuid tensor density 
is isotropic for the fee lattice. 

IV. BCC LATTICE 

Part of the phase diagram for sohd ^He, at higher temperatures and near the melting 
line, includes a bcc phase. Given that the transition temperature for superflow in solids is 
well below 0.5 K, whereas the bcc phase does not occur until temperatures near 1.5 K, it is 
unhkely that the bcc phase will exhibit superflow. Nevertheless, calculation of fg should be 
physically revealing. We anticipate that for a given a/d, since this much more open lattice 
has a lower density and more space, on average, between sites, that the superfiuid fraction 
typically will be lower than for the corresponding fee and hep cases. Calculations on fg for 
this lattice, which has an isotropic superfiuid density, bear this out. 

A. Lattice and Results 

For the bcc lattice, the real-space basis is taken to be 

ai = (a/2)(-l,l,l), a2 = (a/2)(l,-l,l), 03 = (a/2)(l, 1, -1). (17) 

— * 

Then R — iai+ ja2 + ka-^ gives the lattice sites. The reciprocal space basis set is given by 

61 = (27r/a)(0,l,l), 62 = (27r/a)(l,0,l), 63 = (27r/a)(l, 1, 0). (18) 

Then the RLVs are G — mbi + 7162 + n-bs. The bcc unit volume is T^^^cc — (J'tcJ'^i 
Po,bcc = '^/(4cc- Moreover, d = so po,6cc = 3\/3/Ad^ = 1.299/d^. Both the fee and 

hep lattices have the higher densities Po,/cc = Pomp ~ 

V2/d^ = 1.4142/(^3. rp^^ig jjj ^Yiows 
some of our results that are to be compared with Table I for the hep and fee cases. Note 
the convergence rate and note that for the bcc lattice to have the same value of fg as for 
the hep and fee lattices, a/d must be larger. 

Figure 1 shows the superfiuid fraction vs a/d for the bcc lattice. It is typically less than 
that for the fee and bcc lattices, as expected. 
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TABLE III: Superfluid fraction for bcc lattice with a/d = 0.35 {h/d = 0.4950) and two flow 
directions. 



N 


/.(lOO) 


/s(lll) 


2 


0.90084901 


0.90059681 


3 


0.90043751 


0.90041712 


4 


0.90040634 


0.90040469 


5 


0.90040389 


0.90040375 


6 


0.90040369 


0.90040368 



V. CONCLUSIONS 

Given that Kim and Chan obtain SFF's in the range of 1-2%, it would appear that the 
value of h/d is about 0.12.— However, this value is just in the range where the curve of SFF 
vs h/d develops a very steep slope. If a microscopic calculation of h/d (or, better, of pif)) is 
accurate to only 10%, then the actual value of h/d could be overestimated, giving a value for 
fs perhaps as large as 20%. If this is the case, the difference between theory and experiment 
could possibly be attributed to the error in the microscopic calculations of pif), rather than 
to an inadequacy of the one-body theory of superflow. However, current microscopic p{rys 
are expected to be accurate to a few percent, which should enable a judgement to be made 
on whether the one-body phase gives an accurate description of the superfluid. Work is in 
progress to go beyond a one-body phase.-^ 

Finally, we note that, although fg < 1, even at temperature T = 0, this does not imply 
that thermal excitations take up the slack, since we expect that the normal fluid fraction 
/„ = at T = 0. Rather, a lattice fraction must be added to Andreev and Lifshitz's 
macroscopic theory of superflow.— 

VI. ACKNOWLEDGEMENTS 

WMS would like to acknowledge conversations and correspondence with D. Galli and M. 
Chan. This work was supported in part by the Department of Energy, through DOE Grant 
No. DE-FG03-96ER45598. 



* Electronic address: IwsaslowOtamu.ecJul 

1 A. F. Andreev and I. M. Lifshitz, Sov. Phys. JETP 29, 1107 (1969). 

8 



2 A. J. Leggett, Phys. Rev. Lett. 25, 1543 (1970). 

^ A. J. Leggett, Physica Fennica 8, 125 (1973). 

^ A. J. Leggett, J. Stat. Phys. 93, 927 (1998). 

^ W. M. Saslow, Phys. Rev. Lett. 36, 1151 (1976). 

6 E. Kim and M. H. W. Chan, Nature 427, 225 (2004). 

^ E. Kim and M. H. W. Chan, J. Low Temp. Phys. 138, 859 (2005). 

8 E. Kim and M. H. W. Chan, Science 305, 1941 (2004). 

^ E. Kim and M. H. W. Chan, talk at the March 2005 Meeting of the American Physical Society. 

Reported in Science 308, 190 (2005). 

° W. M. Saslow, Phys. Rev. B 71, 092502 (2005). 

^ D. M. Ceperley and B. Bernu, Phys. Rev. Lett. 93, 155303 (2004). 

2 D. E. Galh, M. Rossi, and L. Reatto, Phys. Rev. B 71, 140506 (2005). 

^ D. E. Galli, L. Reatto, and W. M. Saslow, in preparation. 

^ A. J. Leggett, Rev. Mod. Phys. 71, S318 (1999). 

^ C. N. Yang, Rev. Mod. Phys. 34, 644 (1962); F. Bloch and H. E. Rorshach, Phys. Rev.128, 1697 
(1962); W. Kohn, Phys. Rev. 133, A171 (1964). 



9 



